Complete, Theoretical Rovibronic Spectral Characterization of the Carbon Monoxide, Water, and Formaldehyde Cations.

New high-level ab initio quartic force field (QFF) methods are explored which provide spectroscopic data for the electronically excited states of the carbon monoxide, water, and formaldehyde cations, sentinel species for expanded, recent cometary spectral analysis. QFFs based on equation-of-motion ionization potential (EOM-IP) with a complete basis set extrapolation and core correlation corrections provide assignment for the fundamental vibrational frequencies of the A˜2B1 and B˜2A1 states of the formaldehyde cation; only three of these frequencies have experimental assignment available. Rotational constants corresponding to these vibrational excitations are also provided for the first time for all electronically excited states of both of these molecules. EOM-IP-CCSDT/CcC computations support tentative re-assignment of the ν1 and ν3 frequencies of the B˜2B2 state of the water cation to approximately 2409.3 cm-1 and 1785.7 cm-1, respectively, due to significant disagreement between experimental assignment and all levels of theory computed herein, as well as work by previous authors. The EOM-IP-CCSDT/CcC QFF achieves agreement to within 12 cm-1 for the fundamental vibrational frequencies of the electronic ground state of the water cation compared to experimental values and to the high-level theoretical benchmarks for variationally-accessible states. Less costly EOM-IP based approaches are also explored using approximate triples coupled cluster methods, as well as electronically excited state QFFs based on EOM-CC3 and the previous (T)+EOM approach. The novel data, including vibrationally corrected rotational constants for all states studied herein, provided by these computations should be useful in clarifying comet evolution or other remote sensing applications in addition to fundamental spectroscopy.

Ames-2021 CO2 Dipole Moment Surface and IR Line Lists: Toward 0.1% Uncertainty for CO2 IR Intensities.

A highly accurate CO2 ab initio dipole moment surface (DMS), Ames-2021, is reported along with 12C16O2 infrared (IR) intensity comparisons approaching a 1-4‰ level of agreement and uncertainty. The Ames-2021 DMS was accurately fitted from CCSD(T) finite-field dipoles computed with the aug-cc-pVXZ (X = T, Q, 5) basis for C atom and the d-aug-cc-pVXZ (X = T, Q, 5) basis for O atoms, and extrapolated to the one particle basis set limit. Fitting σrms is 3.8 × 10-7 au for 4443 geometries below 15 000 cm-1. The corresponding IR intensity, SAmes-2021, are computed using the Ames-2 potential energy surface (PES), which is the best PES available for CO2. Compared to high accuracy IR studies for 2001i-00001 and 3001i-00001 bands, SAmes-2021 matches NIST experiment-based intensities [SNIST-HIT16 or SHIT20] to -1.0 ± 1.3‰, or matches DLR experiment-based intensities [SDLR-HIT16/UCL/Ames] to 1.9 ± 3.7‰. This indicates the systematic deviations and uncertainties have been significantly reduced in SAmes-2021. The SUCL2015 (or SHITRAN2016) have larger deviations (vs SDLR) and uncertainties (vs SDLR, SNIST) which are attributed to the less accurate Ames-1 PES adopted in UCL-296 line list calculation. The SAmes-2021 intensity of 12C16O2 and 13C16O2 is utilized to derive new absolute 13C/12C ratios for Vienna PeeDee Belemnite (VPDB) with uncertainty reduced by 1/3 or 2/3. Further evaluation of SAmes-2021 intensities are carried out on those CO2 bands discussed in the HITRAN2020 update paper. Consistent improvements and better accuracies are found in band-by-band analysis, except for those bands strongly affected by Coriolis couplings, or very weak bands measured with relatively larger experimental uncertainties. The Ames-2021 296 K IR line lists are generated for 13 CO2 isotopologues, with 18 000 cm-1 and S296 K > 1 × 10-31 cm/molecule cutoff and then combined with CDSD line positions (except 14C16O2). The Ames-2021 DMS and 296 K IR line lists represent a major improvement over previous CO2 theoretical IR intensity studies, including Ames-2016, UCL-296, and recent UCL DMS 2021 update. A real 1 permille level of agreement and uncertainty will definitely require both more accurate PES and more accurate DMS.

What It Takes to Compute Highly Accurate Rovibrational Line Lists for Use in Astrochemistry.

ConspectusWe review the Best Theory + Reliable High-Resolution Experiment (BTRHE) strategy for obtaining highly accurate molecular rovibrational line lists with InfraRed (IR) intensities. The need for highly accurate molecular rovibrational line lists is twofold: (a) assignment of the many rovibrational lines for common stable molecules especially those that exhibit a large amplitude motion, such as NH3, or have a high density of states such as SO2; and (b) characterization of the atmospheres of exoplanets, which will be one of the main areas of research in astronomy in the coming decades. The first motivation arises due to the need to eliminate lines due to common molecules in an astronomical observation in order to identify lines from new molecules, while the second motivation arises due to the need to obtain accurate molecular opacities in order to characterize the atmosphere of an exoplanet. The BTRHE strategy first consists of using high-quality ab initio quantum-chemical methods to obtain a global potential energy surface (PES) and dipole moment surface (DMS) that contains the proper physics. The global PES is then refined using a subset of the reliable high-resolution experimental data. The refined PES then gives energy-level predictions to an accuracy similar to the reproduction accuracy of the experimental data used in the refinement step in the interpolation region (i.e., within the range of the experimental data used in the refinement step). The accuracy of the energy levels will slowly degrade as they are extrapolated to spectral regions beyond the high-resolution experimental data used in the refinement step. However, because the degradation is slow, the predicted energy levels can be used to assign new high-resolution experiments, and the data from these can then be used in a subsequent refinement step. In this way, the global PES eventually can yield highly accurate energy levels for all desired spectral regions including to very high energies and high J values. We show that IR intensities computed with the BTRHE rovibrational wave functions and the DMS can be very accurate provided one has minimized the fitting error of the DMS and tested the completeness of the DMS. Some examples of our work on NH3, CO2, and SO2 are given to highlight the usefulness of the BTRHE strategy and to provide ideas on how to further improve its predictive power in the future. In particular, it is shown how successive refinement steps, once new high-resolution data are available, can lead to PESs that yield highly accurate transition energies to larger spectral regions. The importance of including nonadiabatic corrections to reduce the J-dependence of errors for H-containing molecules is shown with work on NH3. Another very important aspect of the BTRHE approach is the consistency across isotopologues, which allows for highly accurate line lists for any isotopologue once one is obtained for the main isotopologue (which has more high-resolution data available for refinement).

Highly Accurate Quartic Force Field and Rovibrational Spectroscopic Constants for the Azirinyl Cation (c-C2NH2+) and Its Isomers.

The azirinyl cation is an aromatic cyclic molecule that is isoelectronic with cyclopropenylidene, c-C3H2, and c-C3H3+. Cyclopropenylidene has been shown to be ubiquitous, existing in many different astrophysical environments. Given the similar chemistry between C and N, and the relative abundances between C and N in astrophysical environments, it is expected that there should be aromatic ringed molecules that incorporate N in the ring, but as yet, no such molecule has been identified. To address this issue, the present study uses high levels of electronic structure theory to compute a highly accurate quartic force field (QFF) for the azirinyl cation and its two lowest lying isomers, the cyanomethyl and isocyanomethyl cations. The theoretical approach uses the singles and doubles coupled-cluster method that includes a perturbative correction for connected triple excitations, CCSD(T), together with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core-correlation. The QFF is then used in a second-order vibrational perturbation theory analysis (VPT2) to compute the fundamental vibrational frequencies and rovibrational spectroscopic constants for all three C2NH2+ isomers. The reliability of the VPT2 vibrational frequencies is tested by comparison to vibrational configuration interaction (VCI) calculations, and excellent agreement is found between the two approaches. Fundamental vibrational frequencies and rovibrational spectroscopic constants for all singly substituted 13C, 15N, and D isotopologues are also reported. It is expected that the highly accurate spectroscopic data reported herein will be useful in the identification of these cations in high-resolution experimental or astronomical observations.

The anharmonic quartic force field infrared spectra of hydrogenated and methylated PAHs.

Polycyclic aromatic hydrocarbons (PAHs) have been shown to be ubiquitous in a large variety of distinct astrophysical environments and are therefore of great interest to astronomers. The majority of these findings are based on theoretically predicted spectra, which make use of scaled DFT harmonic frequencies for band positions and the double harmonic approximation for intensities. However, these approximations have been shown to fail at predicting high-resolution gas-phase infrared spectra accurately, especially in the CH-stretching region (2950-3150 cm-1, 3 µm). This is particularly worrying for the subset of hydrogenated or methylated PAHs to which astronomers attribute the observed non-aromatic features that appear in the CH-stretching region of spectral observations of the interstellar medium (ISM). In our previous work, we presented the anharmonic theoretical spectra of three linear PAHs and five non-linear PAHs, demonstrating the importance of including anharmonicities into theoretical calculations. In this work we extend these techniques to two methylated PAHs (9-methylanthracene, and 9,10-dimethylanthracene) and four hydrogenated PAHs (9,10-dihydroanthracene, 9,10-dihydrophenanthrene, 1,2,3,4-tetrahydronaphthalene, and 1,2,3,6,7,8-hexahydropyrene) in order to better understand the aliphatic IR features of substituted PAHs. The theoretical spectra are compared with the spectra obtained under matrix isolation low-temperature conditions for the full vibrational fundamental range and under high-resolution, low-temperature gas-phase conditions for the CH-stretching region. Excellent agreement is observed between the theoretical and high-resolution experimental spectra with a deviation of 0.00% ± 0.17%, and changes to the spectra of PAHs upon methylation and hydrogenated are tracked accurately and explained.

Towards completing the cyclopropenylidene cycle: rovibrational analysis of cyclic N3+, CNN, HCNN+, and CNC.

The simple aromatic hydrocarbon, cyclopropenylidene (c-C3H2), is a known, naturally-occurring molecule. The question remains as to whether its isoelectronic, cyclic, fellow aromatics of c-N3+, c-CNN, HCNN+, and c-CNC- are as well. Each of these are exciting objects for observation of Titan, and the rotational constants and vibrational frequencies produced here will allow for remote sensing of Titan's atmosphere or other astrophysical or terrestrial sources. None of these four aromatic species are vibrationally strong absorbers/emitters, but the two ions, HCNN+ and c-CNC-, have dipole moments of greater than 3 D and 1 D, respectively, making them good targets for rotational spectroscopic observation. Each of these molecules is shown here to exhibit its own, unique vibrational properties, but the general trends put the vibrational behavior for corresponding fundamental modes within close ranges of one another, even producing nearly the same heavy atom, symmetric stretching frequencies for HCNN+ and c-C3H2 at 1600 cm-1. The c-N3+ cation is confirmed to be fairly unstable and has almost no intensity in its ν2 fundamental. Hence, it will likely remain difficult to characterize experimentally.

The anharmonic quartic force field infrared spectra of five non-linear polycyclic aromatic hydrocarbons: Benz[a]anthracene, chrysene, phenanthrene, pyrene, and triphenylene.

The study of interstellar polycyclic aromatic hydrocarbons (PAHs) relies heavily on theoretically predicted infrared spectra. Most earlier studies use scaled harmonic frequencies for band positions and the double harmonic approximation for intensities. However, recent high-resolution gas-phase experimental spectroscopic studies have shown that the harmonic approximation is not sufficient to reproduce experimental results. In our previous work, we presented the anharmonic theoretical spectra of three linear PAHs, showing the importance of including anharmonicities into the theoretical calculations. In this paper, we continue this work by extending the study to include five non-linear PAHs (benz[a]anthracene, chrysene, phenanthrene, pyrene, and triphenylene), thereby allowing us to make a full assessment of how edge structure, symmetry, and size influence the effects of anharmonicities. The theoretical anharmonic spectra are compared to spectra obtained under matrix isolation low-temperature conditions, low-resolution, high-temperature gas-phase conditions, and high-resolution, low-temperature gas-phase conditions. Overall, excellent agreement is observed between the theoretical and experimental spectra although the experimental spectra show subtle but significant differences.

The anharmonic quartic force field infrared spectra of three polycyclic aromatic hydrocarbons: Naphthalene, anthracene, and tetracene.

Current efforts to characterize and study interstellar polycyclic aromatic hydrocarbons (PAHs) rely heavily on theoretically predicted infrared (IR) spectra. Generally, such studies use the scaled harmonic frequencies for band positions and double harmonic approximation for intensities of species, and then compare these calculated spectra with experimental spectra obtained under matrix isolation conditions. High-resolution gas-phase experimental spectroscopic studies have recently revealed that the double harmonic approximation is not sufficient for reliable spectra prediction. In this paper, we present the anharmonic theoretical spectra of three PAHs: naphthalene, anthracene, and tetracene, computed with a locally modified version of the SPECTRO program using Cartesian derivatives transformed from Gaussian 09 normal coordinate force constants. Proper treatments of Fermi resonances lead to an impressive improvement on the agreement between the observed and theoretical spectra, especially in the C-H stretching region. All major IR absorption features in the full-scale matrix-isolated spectra, the high-temperature gas-phase spectra, and the most recent high-resolution gas-phase spectra obtained under supersonically cooled molecular beam conditions in the CH-stretching region are assigned.

Linear transformation of anharmonic molecular force constants between normal and Cartesian coordinates.

A full derivation of the analytic transformation of the quadratic, cubic, and quartic force constants from normal coordinates to Cartesian coordinates is given. Previous attempts at this transformation have resulted in non-linear transformations; however, for the first time, a simple linear transformation is presented here. Two different approaches have been formulated and implemented, one of which does not require prior knowledge of the translation-rotation eigenvectors from diagonalization of the Hessian matrix. The validity of this method is tested using two molecules H2O and c-C3H2D(+).

Quartic force field rovibrational analysis of protonated acetylene, C2H3(+), and its isotopologues.

Protonated acetylene, C2H3(+), is among the simplest carbocations. Comprehensive experimental or highly accurate computational spectroscopic data is lacking for this system due to its inherent complexities. Utilizing state-of-the-art quartic force fields (QFFs), the spectroscopic constants and fundamental vibrational frequencies are provided in this work for the nonclassical, bridged, cyclic global minimum. The rotational constants match experiment to better than 0.1%, and the computed ν2 antisymmetric HCCH stretch is less than 3.0 cm(­1) different from experiment. Hence, the rovibrational spectroscopic data provided herein for c-C2H3(+) and its deuterated isotopologues enrich the chemical understanding of this system. Unfortunately, the same rovibrational spectroscopic data is not as trustworthy for the classical, linear form of protonated acetylene due to the shallow well in which it resides on the potential energy surface. However, spectroscopic data are provided for this isomer in the Supporting Information to enhance future studies.

Highly accurate potential energy surface, dipole moment surface, rovibrational energy levels, and infrared line list for ³²S¹6O2 up to 8000 cm⁻¹.

A purely ab initio potential energy surface (PES) was refined with selected (32)S(16)O2 HITRAN data. Compared to HITRAN, the root-mean-squares error (σ(RMS)) for all J = 0-80 rovibrational energy levels computed on the refined PES (denoted Ames-1) is 0.013 cm(-1). Combined with a CCSD(T)/aug-cc-pV(Q+d)Z dipole moment surface (DMS), an infrared (IR) line list (denoted Ames-296K) has been computed at 296 K and covers up to 8000 cm(-1). Compared to the HITRAN and CDMS databases, the intensity agreement for most vibrational bands is better than 85%-90%. Our predictions for (34)S(16)O2 band origins, higher energy (32)S(16)O2 band origins and missing (32)S(16)O2 IR bands have been verified by most recent experiments and available HITRAN data. We conclude that the Ames-1 PES is able to predict (32/34)S(16)O2 band origins below 5500 cm(-1) with 0.01-0.03 cm(-1) uncertainties, and the Ames-296K line list provides continuous, reliable and accurate IR simulations. The K(a)-dependence of both line position and line intensity errors is discussed. The line list will greatly facilitate SO2 IR spectral experimental analysis, as well as elimination of SO2 lines in high-resolution astronomical observations.

Fundamental vibrational frequencies and spectroscopic constants of cis- and trans-HOCS, HSCO, and isotopologues via quartic force fields.

Highly accurate, coupled-cluster-based quartic force fields (QFFs) have been employed recently to provide spectroscopic reference for a myriad of molecules. Here, we are extending the same approach to provide vibrational and rotational spectroscopic reference data for the sulfur analogues of HOCO, HSCO, and HOCS, in both the cis and trans conformations as well as the D and (34)S isotopologues of each system. The resulting energies corroborate previous computations showing that trans-HSCO is the lowest-energy isomer for this system. The vibrational frequencies are computed with both second-order vibrational perturbation theory (VPT2) and vibrational configuration interaction (VCI) methods. The VPT2 and VCI QFF frequencies largely agree with one another to better than 5.0 cm(-1) (often better than 1.0 cm(-1)) and are also consistent with the type of behavior exhibited in previous studies. As such, the reference data provided here should assist in analysis of environments in which these sulfur systems may be found, including the interstellar medium, combustion flames, or laboratory simulations of either.

Limited rotational and rovibrational line lists computed with highly accurate quartic force fields and ab initio dipole surfaces.

In this work, computational procedures are employed to compute the rotational and rovibrational spectra and line lists for H2O, CO2, and SO2. Building on the established use of quartic force fields, MP2 and CCSD(T) Dipole Moment Surfaces (DMSs) are computed for each system of study in order to produce line intensities as well as the transition energies. The computed results exhibit a clear correlation to reference data available in the HITRAN database. Additionally, even though CCSD(T) DMSs produce more accurate intensities as compared to experiment, the use of MP2 DMSs results in reliable line lists that are still comparable to experiment. The use of the less computationally costly MP2 method is beneficial in the study of larger systems where use of CCSD(T) would be more costly.

Anharmonic rovibrational calculations of singlet cyclic C4 using a new ab initio potential and a quartic force field.

We report a CCSD(T)/cc-pCV5Z quartic force field (QFF) and a semi-global CCSD(T)-F12b/aug-cc-pVTZ potential energy surface (PES) for singlet, cyclic C4. Vibrational fundamentals, combinations, and overtones are obtained using vibrational second-order perturbation theory (VPT2) and the vibrational configuration-interaction (VCI) approach. Agreement is within 10 cm(-1) between the VCI calculated fundamentals on the QFF and PES using the MULTIMODE (MM) program, and VPT2 and VCI results agree for the fundamentals. The agreement between VPT2-QFF and MM-QFF results is also good for the C4 combinations and overtones. The J = 1 and J = 2 rovibrational energies are reported from both VCI (MM) on the PES and VPT2 on the QFF calculations. The spectroscopic constants of (12)C4 and two C2v-symmetry, single (13)C-substituted isotopologues are presented, which may help identification of cyclic C4 in future experimental analyses or astronomical observations.

Protonated nitrous oxide, NNOH+: fundamental vibrational frequencies and spectroscopic constants from quartic force fields.

The interstellar presence of protonated nitrous oxide has been suspected for some time. Using established high-accuracy quantum chemical techniques, spectroscopic constants and fundamental vibrational frequencies are provided for the lower energy O-protonated isomer of this cation and its deuterated isotopologue. The vibrationally-averaged B0 and C0 rotational constants are within 6 MHz of their experimental values and the D(J) quartic distortion constants agree with experiment to within 3%. The known gas phase O-H stretch of NNOH(+) is 3330.91 cm(-1), and the vibrational configuration interaction computed result is 3330.9 cm(-1). Other spectroscopic constants are also provided, as are the rest of the fundamental vibrational frequencies for NNOH(+) and its deuterated isotopologue. This high-accuracy data should serve to better inform future observational or experimental studies of the rovibrational bands of protonated nitrous oxide in the interstellar medium and the laboratory.

The 1 3A' HCN and 1 3A' HCO+ vibrational frequencies and spectroscopic constants from quartic force fields.

Building on previous studies involving coupled cluster quartic force fields for the description of spectroscopic constants and vibrational frequencies of astronomically relevant molecules, this work applies the same techniques to the elucidation of such properties for the bent 1 (3)A' state of HCN and the isoelectronic 1 (3)A' HCO(+). Core correlation is treated both by explicit means and as a correction. Each approach gives closely comparable spectroscopic constants and vibrational frequencies once more, indicating that the composite method is a viable and less costly alternative. We are providing fundamental vibrational frequencies for these systems where agreement with experiment in previous studies has been within 4 cm(-1) or better. Frequencies for the first overtones and combination bands as well as various spectroscopic constants are also reported.

Dipole surface and infrared intensities for the cis- and trans-HOCO and DOCO radicals.

The vibrational spectra for the HOCO radical in both of its conformers and deuterated isotopologues are shown here for the first time. Building on previous work with coupled cluster quartic force fields (QFFs) in the computation of the fundamental vibrational frequencies for both cis- and trans-HOCO, coupled cluster dipole surfaces are now provided for both HOCO conformers and their corresponding deuterated isotopologues. These surfaces and subsequent vibrational configuration interaction (VCI) computations produce the intensities of transitions into vibrational states including the fundamentals, overtones, and first few combination bands of less than 4000 cm(-1), slightly beyond the O-H stretch. Simulated spectra with an artificial full width at half-maximum broadening of 10 cm(-1) are also provided in order to aid in the characterization of HOCO's vibrational frequencies and to assist detection in various laboratory or astronomical observations.

Fundamental vibrational frequencies and spectroscopic constants of HOCS+, HSCO+, and isotopologues via quartic force fields.

Besides the ν(1) O-H stretching mode at 3435 cm(-1) for HOCS(+), the fundamental vibrational frequencies for this cation and its HSCO(+) isomer have not been determined experimentally. Because these systems are analogues to HOCO(+), a detected interstellar molecule, and are believed to play an important role in reactions of OCS, which has also been detected in the interstellar medium, these cations are of importance to interstellar chemistry and reaction surface studies. This work provides the fundamental vibrational frequencies and spectroscopic constants computed with vibrational perturbation theory (VPT) at second order and the vibrational configuration interaction (VCI) method conjoined with the most accurate quartic force field (QFF) applied to date for these systems. Our computations match experiment to better than 2 cm(-1) for the known O-H stretch. Additionally, there is strong agreement in the prediction of the fundamentals across methods and choices of QFFs. The consistency in the computations and the correspondence for the known mode should give accurate reference data for the rovibrational spectra of these cations and their singly substituted isotopologues for D, (18)O, and (34)S.

Quartic force field predictions of the fundamental vibrational frequencies and spectroscopic constants of the cations HOCO+ and DOCO+.

Only one fundamental vibrational frequency of protonated carbon dioxide (HOCO(+)) has been experimentally observed in the gas phase: the ν(1) O-H stretch. Utilizing quartic force fields defined from CCSD(T)/aug-cc-pVXZ (X = T,Q,5) complete basis set limit extrapolated energies modified to include corrections for core correlation and scalar relativistic effects coupled to vibrational perturbation theory and vibrational configuration interaction computations, we are predicting the full set of gas phase fundamental vibrational frequencies of HOCO(+). Our prediction of ν(1) is within less than 1 cm(-1) of the experimental value. Our computations also include predictions of the gas phase fundamental vibrational frequencies of the deuterated form of the cation, DOCO(+). Additionally, other spectroscopic constants for both systems are reported as part of this study, and a search for a cis-HOCO(+) minimum found no such stationary point on the potential surface indicating that only the trans isomer is stable.

An isotopic-independent highly accurate potential energy surface for CO2 isotopologues and an initial (12)C(16)O2 infrared line list.

An isotopic-independent, highly accurate potential energy surface (PES) has been determined for CO(2) by refining a purely ab initio PES with selected, purely experimentally determined rovibrational energy levels. The purely ab initio PES is denoted Ames-0, while the refined PES is denoted Ames-1. Detailed tests are performed to demonstrate the spectroscopic accuracy of the Ames-1 PES. It is shown that Ames-1 yields σ(rms) (root-mean-squares error) = 0.0156 cm(-1) for 6873 J = 0-117 (12)C(16)O(2) experimental energy levels, even though less than 500 (12)C(16)O(2) energy levels were included in the refinement procedure. It is also demonstrated that, without any additional refinement, Ames-1 yields very good agreement for isotopologues. Specifically, for the (12)C(16)O(2) and (13)C(16)O(2) isotopologues, spectroscopic constants G(v) computed from Ames-1 are within ±0.01 and 0.02 cm(-1) of reliable experimentally derived values, while for the (16)O(12)C(18)O, (16)O(12)C(17)O, (16)O(13)C(18)O, (16)O(13)C(17)O, (12)C(18)O(2), (17)O(12)C(18)O, (12)C(17)O(2), (13)C(18)O(2), (13)C(17)O(2), (17)O(13)C(18)O, and (14)C(16)O(2) isotopologues, the differences are between ±0.10 and 0.15 cm(-1). To our knowledge, this is the first time a polyatomic PES has been refined using such high J values, and this has led to new challenges in the refinement procedure. An initial high quality, purely ab initio dipole moment surface (DMS) is constructed and used to generate a 296 K line list. For most bands, experimental IR intensities are well reproduced for (12)C(16)O(2) using Ames-1 and the DMS. For more than 80% of the bands, the experimental intensities are reproduced with σ(rms)(ΔI) < 20% or σ(rms)(ΔI∕δ(obs)) < 5. A few exceptions are analyzed and discussed. Directions for future improvements are discussed, though it is concluded that the current Ames-1 and the DMS should be useful in analyzing and assigning high-resolution laboratory or astronomical spectra.